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Abstract. Transport properties of degenerate relativistic 
electrons along quantizing magnetic fields in neutron star 
crusts are considered. A kinetic equation is derived for the 
spin polarization density matrix of electrons. Its solution 
does not depend on the choice of basic electron wave func- 
tions unlike previous solutions of the traditional kinetic 
equation for the distribution function. The density matrix 
formalism shows that one can always reach high accuracy 
with the traditional method by a proper choice of the basic 
functions. Electron Coulomb scattering on ions is consid- 
ered in liquid matter, and on high-temperature phonons or 
on charged impurities in solid matter. In the solid regime, 
the Debye - Waller reduction of phonon scattering can 
strongly enhance the longitudinal thermal or electric con- 
ductivity. An efficient numerical method is proposed for 
calculating the transport properties of electron gas at any 
magnetic field of practical interest. 

Key words: stars: neutron - dense matter - magnetic 
fields - radiation mechanisms: thermal 



1. Introduction 

Accurate transport coefficients in neutron star crusts are 
important for analysing the thermal evolution of neutron 
stars and evolution of their magnetic fields. In outer crusts 
of cooling magnetized neutron stars, the heat is mainly 
transported along the magnetic fields. There exist several 
competing heat transport mechanisms across the field, but 
the longitudinal currents are carried mostly by electrons 
through their scattering on phonons or charged impurities 
in the solid phase and on ions in the liquid phase. As a 
rule, the electrons in the crust are strongly degenerate 
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and may be relativistic; the magnetic held can be easily 
quantizing. Transport properties of the crusts have been 
studied in a number of papers (e.g., Yakovlev & Kaminker 
1994, and references therein). 

In the present work, the most important problem of 
longitudinal electron transport in quantizing magnetic 
helds is studied with the use of the quantum density ma- 
trix formalism, instead of the traditional kinetic equation 
for the electron distribution function employed in previ- 
ous studies. The main advantage of the present approach 
is that it is independent of the choice of the basis of elec- 
tron states (basis states are not unique due to the electron 
spin degeneracy). 

We consider three main electron scattering mecha- 
nisms. The first one is the Coulomb scattering on ions 
in the liquid or gaseous phase. The second one is the scat- 
tering on high-temperature phonons in the solid phase. In 
the latter case, we take into account the Debye - Waller 
factor whose importance has been emphasized and proved 
by Itoh et al. (1984b, 1993) for the non-magnetic case. 
We show that the effect of this factor is much stronger in 
quantizing magnetic fields. The third mechanism is the 
Coulomb scattering on charged impurities in the solid 
phase, important much below the melting temperature. 

The paper is composed as follows. In Sect. || we de- 
scribe the physical conditions of interest, electron scat- 
tering potentials, transport coefficients and their expres- 
sions in the non-magnetic case. In Sect. [}] we derive a lin- 
earized kinetic equation for the density matrix and com- 
pare its numerical solutions with the traditional solutions 
employed in all previous works. Mathematical properties 
of the new equation are discussed in Appendix A. The ef- 
fect of the Debye - Waller factor in quantizing magnetic 
helds is studied in Sect. [|. The results are summarized 
in Sect. ^[ In Appendix B we present new expressions for 
some intermediate integrals. These expressions ensure effi- 
cient computation of the transport properties for the case 
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when many Landau orbitals are occupied. Previous results 
(Yakovlev 1984, Hernquist 1984, Schaaf 1988, Van Riper 
1988) were restricted to 30 Landau orbitals at most. 

2. Basic equations 

2.1. Physical conditions 

Magnetic fields B in magnetized neutron star crusts are 
known to range between 10 11 and 10 14 G. For tempera- 
tures T ~ (10 6 - 10 s ) K and densities < 10 8 g cm -3 of 
interest, the Landau quantization of electron motion in 
a magnetic field can be very important (then the field is 
called quantizing). Consider not too low densities (typi- 
cally, ^ 10 5 g cm~ 3 at B ~ 10 12 G, see, e.g., Yakovlev 
1984, Van Riper 1988) when the electron gas is almost 
free and the ions are fully ionized. An electron energy is 
then given by 



1 /2 

i(Pz) = c ((m c c) 2 + 2Uuj B m c n+p 2 z ) 



(1) 



where p z is the electron momentum along the field, ojb = 
eB/(m c c) is the electron cyclotron frequency, (— e) is the 
electron charge, m e is the electron mass, c is the speed of 
light, and n = 0, 1, 2, ... is a Landau quantum number. 
The electron number density can be expressed as 



m c u>B 
(2nh) 2 
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(2) 



where s is a spin quantum number (s = —1 for n = 0, and 
s = ±1 for n = 1, 2, . . ., see Sect. 0), 



Me) = 



exp 



(3) 



is the Fermi - Dirac distribution function, and \i is the 
chemical potential (including m c c 2 ). A fitting formula for 
calculation of n c is given in Appendix C. 

The state of ions is mainly determined by the ion cou- 
pling parameter T = (Z "e) 2 / (ak^T) , where Ze is an ion 
charge, /cb is the Boltzmann constant, a — (47rni/3)~ 1//3 
is the ion sphere radius, and n\ = n c /Z is the ion num- 
ber density. If T <C 1, the ions constitute an ideal gas. 
For higher T, the gas gradually transforms into a strongly 
coupled Coulomb liquid. The liquid solidifies at the melt- 
ing temperature T = T m which corresponds to T = 172 
(Nagara et al. 1987). The zero-point quantum vibrations 
of the Coulomb crystal become comparable with the ther- 
mal vibrations at k^T hu> p , where w p = \J ' ArKZ 2 e 2 n\jm\ 
is the ion plasma frequency, and m\ is the ion mass. We 
will not consider superstrong magnetic fields (discussed, 
e.g., by Yakovlev 1984) which affect the properties of the 
ion plasma component. 

A more detailed description of the physical conditions 
is given by Yakovlev (1984) and Van Riper (1988). 



2.2. Scattering potentials 

Consider three important cases when the electron scatter- 
ing is almost elastic. The first case is the Coulomb scatter- 
ing on ions in the liquid or gaseous phase (T > T m ). The 
second case is the scattering on high-temperature phonons 
(TiWp/fcB < T m ). The third one is the Coulomb scat- 
tering on charged impurities in the lattice, which is im- 
portant for k^T <C hujp. The impurities represent ions of 
charge Z- lmp ^= Z immersed accidentally in lattice sites. 
Thus our results will cover a wide range of temperatures. 

The electron- ion scattering potential V(r) in the liq- 
uid or gaseous phase can approximately be taken as a 
screened Coulomb potential (Yakovlev 1984). Its Fourier 
image U(q) is given by 



[4nZe 2 /( q 2 + r; 2 )] 2 . 



(4) 



where r s is an effective screening length, r~ 2 = + 
r~ 2 . Here r; and r e are the screening lengths due to ions 
and electrons, respectively. In the most important liquid 
regime (1 ^ T < 172), according to Yakovlev (1984), the 
ion screening length is 



(5) 



fi = aye/6 « 0.67a 



(Hernquist (1984) used a less accurate approximation for 
r;). The electron screening length is determined as 



[47te 2 dn e /dfi 



-1/2 



(G) 



When the temperature is low enough (fceT <ti Tioj* b with 
u* B = UB/n e c 2 //i), we have 



(a m /n) 2 

(a m /r c ) 2 



^ {m c hLo B ) 



-1/3 



3ttZ 



-i 2/3 



2a 



E 



\Pz\c 



(7) 
(8) 



where a m = (hc/eB) 1 / 2 is the magnetic length, a = e 2 /hc 
is the fine-structure constant, and the energy and momen- 
tum variables are assumed to be taken on the Fermi sur- 
face: e = t n {p z ) — (i. At arbitrary temperature, the fitting 
formula of Appendix C can be used to calculate n c and 
dn e /dfi for estimation of r; and r c according to Eqs. (|^) 
and (|). 

For scattering on high-temperature phonons in the 
solid phase, one has 



IMq)|' 
where 



exp [-2W(q)} 



3fc B Tu_ 2 
4:Tm\Z 2 e 2 



u- 2 a 



(9) 



(10) 
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is the mean squared thermal displacement of ions, u-2 
is a numerical factor determined by the phonon spec- 
trum (it_2 = 13 for the bcc lattice), and e~ 2W is the 
Debye - Waller factor. The latter factor is usually neg- 
ligible for scattering in terrestrial solids (e.g., Davydov 
1976), but it is important in dense neutron star matter 
(Itoh et al. 1984b, 1993). For the high-temperature solids 
(k^T J> TiLO-p) of interest, one has (e.g., Itoh et al. 1984a) 



2VF(q) « (r T g) 2 /3. 



(11) 



Note that Eq. (j^) represents the familiar high-temperature 
asymptote of the one-phonon scattering potential (Yakov- 
lev & Urpin 1980) multiplied by the Debye - Waller term 
(Itoh et al. 1984b) to include multiphonon processes. 

Finally, the Coulomb scattering on impurities corre- 
sponds to (e.g., Yakovlev & Urpin 1980) 



|[/ imp (q)| 2 = [MZimp - Z)e 2 /(q 2 + r' 2 )]' 



(12) 



In this case the screening length r s is most likely deter- 
mined by the electrons: r s = r c . The Coulomb scattering 
on impurities is very similar to that on ions, and we do 
not consider them separately in detail. The results for im- 
purities can be obtained from those for ions by replacing 



where a is the longitudinal electric conductivity. 
For practical use, Eq. (17) can be rewritten as 



e dz a a dz 



q=-^Tj-J£, (18) 
a az 



where (3 /a and ac = A — T0 2 /a are the longitudinal ther- 
mopower and thermal conductivity, respectively. 

For nearly elastic electron scattering, the kinetic coef- 
ficients er, fS, and A may be expressed as 




e(e-/i)/T 



n c r(e)c 2 / <9/o 



'if ) de, (19) 



where r(e) is the effective energy-dependent relaxation 
time for the electrons. 

2.4- Non-magnetic electron relaxation times 

In the absence of the magnetic field (or for non-quantizing 
field) , the inverse effective relaxation time (effective colli- 
sion frequency) of an electron with energy e can be pre- 
sented as 



L (e) = 7ii«otr(e), 



(20) 



Z ► Zi mp — Z, n\ 



(13) 



where v is the electron velocity and <x tr (e) is the transport 
cross section: 



where nj mp is the impurity number density. 

Equations (|j) and (^) can be conveniently written as 



I Cic 



, v ,2 2-7T 

i(q) = —; 

nil 



2hc/a B 



|tfph(q)r 



7T (2hc 
ml \ q 



-2W 



where I is a scale length (Yakovlev 1984) 
■" " Tuor , 3 / he 



l\on 



27rniZ 2 e 4: 



47mi \Ze 2 rT 



For the scattering on impurities, Ui mp and I 



nmp 



(14) 
(15) 

(16) 
are ob- 



tained from [/ion and Zi on by using Eq. (|T3|). 
2.3. Transport coefficients 

Let j and q be the densities of the electric and thermal cur- 
rents induced by sufficiently weak electric field £ and gra- 
dients of temperature T and electron chemical potential /j, 
directed along the magnetic field (along the z-axis). The 
currents are determined by three transport coefficients a, 
P and A, 



J = ° 



e oz 



~(3T £ 



1 dfi 

e dz 



A 



dT 



(17) 



47T 



dfi'cr(p^p')(l-COSe). 



(21) 



Here p and p' are electron momenta before and after scat- 
tering, respectively, dSl and dfl' are solid angle elements, 
is the scattering angle, and er(p — > p') is a differential 
scattering cross section. In the Born approximation, 



cr(p 



2^2 



l^(q)l 2 e 

47T 2 ?i 4 C 4 



e 



* Sm 2 



(22) 



Consider the Coulomb scattering. Let us substitute the 
Fourier image (||) into Eq. (|2^) . Integrating in Eq. ( pl| ) , we 
arrive at the well known result (e.g., Yakovlev and Urpin 
1980): 



o-tr(e) = 4tt ( — ) A(c). 
pv 



(23) 



Here A(e) is the Coulomb logarithm (Yakovlev 1980): 
A(e) = \ [ln(l + W )-(l + W - 1 )- 1 ] - 



v 

2c 2 



1 - 2w^ 1 In (l + w) + (l + w)~ 



(24) 



and 



= (2pr s /h) 2 . The terms of order u> -1 are of- 
ten neglected in A(e) (e.g., Yakovlev and Urpin 1980, 
Yakovlev 1984). However these terms are very significant 
in Coulomb liquids of carbon and lighter elements. 
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Analogously, for the scattering on phonons, from 
Eqs. M) and (111]) we obtain 



•wW-T^)''*^"-^"^)- (25) 

where 

Ri = - (1 - e' w ) , Jfe = (i_ e -»(l + «,)), (26) 

and to = (4/3) (prr/fr) 2 - Generally, we have R\ < 1 and 
R-2 < 1 due to the Debye - Waller factor. If w -C 1, 
the Debye - Waller factor is insignificant, i?i = R2 = 1, 
and Eq. (^5|) reproduces the result of Urpin and Yakovlev 
(1980). 

Finally, one obtains 



9 

p v 



47rA(e)Z 2 e 4 ni ' 



3?i 2 t> 



87r4Z 2 e 4 7ii 



iil(e) - ifc(e) ^ 



(27) 
(28) 



for the Coulomb and phonon scattering, respectively. 

3. Density matrix 

3.1. Kinetic equation 

Transport properties of a degenerate relativistic electron 
gas in quantizing magnetic fields were studied by Yakovlev 
(1984), Hernquist (1984), Van Riper (1988), and Schaaf 
(1988) on the basis of the linearized kinetic equations 
for the electron distribution function. However, as noted 
by Yakovlev (1984), a selfconsistent description should be 
based on the quantum density matrix formalism. 

Let us consider a uniform magnetic field B along the 
z-axis and use the Landau gauge of the vector potential: 
A = (—By, 0,0). Quantum states of a free electron in 
the magnetic field form a complete orthogonal basis. The 
basic states can be labelled by the quantum numbers K = 
(px)Pzi n , s), where p x determines the y-coordinate of the 
guiding center, ijb = p x /{m e WB)- An explicit solution of 
the Dirac equation reads (e.g., Sokolov & Ternov 1968) 



exp[i(p x x +p z z)/h] 
^(r) = -. u/2 Xn S {Pz,y-VB), 



(29) 



where L x and L z are the normalization lengths, and Xn 
can be chosen as 



Xn,i(Pz,y) 



Xn-l(Pz,y)- 



-1/2 



^/2E(E + 1) 



-1/2 



y/2E(E+l) 



f(E+l)H n -i(y/a m )\ 


P vn {E)H n -i{y/a m ) 
\ -V2bnH n {y/a m ) ) 

I \ 

(E + l)H n {y/a m ) 

m ) 

\-P ri n(E)H n (y/a m ) J 



(30) 



Here b = Tiojb/ (m c c 2 ), E = e/(m c c 2 ) and P vn (E) 
7, (E 2 — 1 — 2on) 1 / 2 are, respectively, the magnetic field, the 
energy and the longitudinal momentum in the relativistic 
units, r) = signp z , 



7r l/4( 2 n n !)l/i 



(31) 



and i?n(0 is a Hermitc polynomial. 

In the non-relativistic limit, the two lower components 
of the bispinors (30) are negligible and this basis corre- 
sponds to fixed spin projections (sh/2) on the z-axis. How- 
ever one can use another basis: 



Xr. 



Xns COS ( 



SXr, 



, sin < 



(32) 



It is sufficient to assume that < <f> < tt/2; <f> may depend 
on n and should vanish for n = 0. In particular, the choice 



arcsm — 
2 



l-(E 2 -l)-V 2 P nn (E) 



1/2 



(33) 



yields the basis of states with fixed helicity used by Hern- 
quist (1984) and Schaaf (1988). 

Using the distribution function formalism, one gen- 
erally obtains different conductivities with different ba- 
sic functions. Yakovlev (1984) employed both </> = and 
4> = 4> n and interpolated between corresponding results 
assuming that the former choice is appropriate in the non- 
relativistic limit (E ~ 1) while the latter one is adequate 
in the ultrarelativistic case (£> 1). 

The accuracy of the distribution function calculations 
can be verified using the density matrix, 



Pk'k = (tp K '\p\lp K ), 



(34) 



where p is the statistical operator. 

Let us consider elastic collisions and assume that the 
density matrix is diagonal in energy: e' = e. Furthermore, 
we assume the diagonality in momenta. We shall show 
that the latter property is not violated by collisions. The 
diagonalities in e and p z lead to the diagonality in n. The 
diagonality of the density matrix in p x , p z allows us to 
treat the dependence of p K / K on x, z parametrically, with- 
out using the Wigner transformation. Thus, we can write 



Pk'k — fijl'Tl ^pLpx ^p' z Pz Pns' sip^ i % , p x , P z ) : 



(35) 



where 5 P > aPa = {2ttTi/ L a ) 5(p' a — p a ). Analysing the lon- 
gitudinal transport properties, we can assume that p is 
independent of x. Let us introduce the spin-polarization 
density matrix, 



p nns > s {z,e) = J p ns , s (z,p x ,p z ) ^-dp x . 



(36) 
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Then we obtain the kinetic equation with the classical 
left-hand side 



d d d 

dt z dz ' dp z 



dt 



(37) 



where £ is a longitudinal electric field and v z is the ve- 
locity. The density matrix depends on p z through e and 
7/. The right-hand side of Eq. (|37]) is the collision integral 
to be determined from microscopic considerations. Owing 
to the linearity of the equations which govern the electron 
wave- function evolution during a scattering event, the col- 
lision integral should be linear in p rinslS2 . Thus it can be 
presented as 



dt 



(38) 



The equilibrium density matrix is /c4il 2 ( e ) = fo( e )^s 1 s 2 , 
and the collision integral vanishes at the equilibrium. 

Deviations from the equilibrium in the linear regime 
can be treated as small perturbations. In the zero-order 
approximation, the density matrix is equal to p^ ; it de- 
pends on z parametrically through /i and T. In the first- 
order approximation, it is customary to write 



Pr\ns\s<i {z, c) 



Pi%(e)+Vl 
e-p. dT 



dfo(e) 



e£ 



dp 

dz 



T dz 



l Prjns 1 s 2 ( e )) 



(39) 



where I is the scale length defined by Eq. ([16]). Keeping the 
zero-order terms on the left-hand side of Eq. ( |37j ) we ob- 
tain the set of algebraic equations for the non-equilibrium 
corrections ip: 



I 



^ ' V -^r;'Ti's' 1 s^;r;nsiS2 Vf] 1 s i s 2 V ^SiS 



(40) 



■q'n' s-^s^ 



A solution of these equations allows us to calculate the 
kinetic coefficients cr, 0, and A in the transport relations 
(17) and so determine the longitudinal electric and ther- 
mal conductivities and thermopower. 

Using Eq. (39) we can write cr, f3 and A in the form of 
Eq. (|ll]) , with the effective relaxation time defined as 



. _ dm c uj B 
2( / nTic) 2 n c 

*W = 2 E Vnnss(e). 

r;=±l n=0 s=±l 



(41) 
(42) 



Here n e = Int(i/) specifies the highest Landau level popu- 
lated by electrons with energy e, and 



e — m r c 2 e + m c c 2 



hu>B 2m e c 2 



(E 2 -l)/(2b). 



(43) 



Equations ( |19[ ) and (41) reproduce Eq. (24) of Yakovlev 
(1984). However now Eq. (42) contains diagonal elements 
of the density matrix (to be determined from Eq. (Eo[)) 
instead of the distribution function in the traditional ap- 
proach. The factors D in Eq. (^fj) are obtained in Sect. 
|3.2| . They are found to be much more complicated than 
analogous factors in the traditional equations. 

3.2. Collision integral 

First let us derive the collision integral for the full density 
matrix p K ' K . After that the collision integral ( |3^ ) is ob- 
tained by summing over those quantum numbers in which 
p K / K is diagonal. 

In the "quasiclassical approach" for the density ma- 
trix (Sobelman et al. 1981), interaction of electrons with 
perturbers (ions or phonons) is described by a scattering 
potential Vo(t). A wave function of an interacting electron 
evolves as tp(t) — S(t,to)tp(to), where S is the scatter- 
ing operator. Substituting this into Eq. ([m|) and averaging 
over collision parameters we obtain 



(44) 



where S K i K = (K,'\S(t, —oo)\k) is the scattering matrix, 
and the brackets (. . .) pc denote space averaging over per- 
turber centers ro (i.e., space integration with the weight 
rij), while the brackets without subscript denote the 
quantum-mechanical averaging. The sum over k' includes 
that over p' x and p' z , which should be performed according 
to the correspondence rule 



E 



L, 



dp x , z - 



(45) 



Following the "adiabatic switch on" method (Landau & 
Lifshitz 1976), we put 



V (t) = V(r-r )e e 



(46) 



where e — » is the adiabatic parameter, and V(r) is the 
actual scattering potential. Now we obtain 



dp K 



dt 
where 



^ ] D k,' 1 k' 2 k, 1 k 2 Pk'jK^i 



(47) 



lim . 

£^o \ dt 



t=o 



(48) 



The scattering matrix elements S K i K can be calculated in 
the second-order perturbation theory (e.g., Landau & Lif- 
shitz 1976). Let us substitute them into Eqs. @ and @. 



(i 
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Then the terms linear in V are canceled after averaging 
over ro, while the quadratic terms give 



D 



K , K. g , K l K 2 



71 [ S ^2 S i u '^ 1 )(V* Kl V KK ' i ) l 



1 

id 



fc2 H 



1 



(49) 



where u) K > K = (e' — e)/h~, and V" K ' K is the matrix element of 
the potential V. In deriving Eq. (49) we have taken into 
account the well known relationship 



lim 



1 



1 



: tt6{uj)±iP-, 
£^0 e =p iui UJ 



(50) 



where V denotes the Cauchy main part. An averaged bi- 
nary product in Eq. (49) is equal to 

27m; f 
(K*i« 1 K^ 2 )pc = , L L \ 2 S(qi x - q2x)S{qiz - q2z)J dq y x 

|[/(q)| 2 M s *, Si {n',p' z ; m,pw, u)M s > 2S2 (n',p' z ;n2,P2 Z ]U),(51) 

where q = (p' p)/h~, U(q) is the Fourier image of the 
potential, and 

M s , s (n',p' z ;n,p z ;u) = 

xt>M,V - q x </2) X ns(Pz,y + q x a 2 j2yiyy&y (52) 



is the matrix element which actually depends on q x and 
q y only through the variable u = (q% + g^)a^/2 (e.g., 
Kaminkcr & Yakovlev 1981). In Eq. (51) we have set 

Pix = P2x> Piz = P2z> and n i = n 2, since the density 
matrix is diagonal in these arguments. Due to the same 
diagonality, u K / Kl — w K ^ K2 , and the last term in Eq. (49) 
vanishes. 

Finally, let us sum D K ^ Kl ^ 2 over n 2) Pi X) P2x, and 
P2 Z using Eq. (^5j) . Thus we arrive at the collision integral 
@, where 



^rj'n'si s' 2 :7]nsiS2 



\P'r}ns\S2 ,7)' n' s^s 2 



Explicit expressions for the real coefficients a and complex 
coefficients A depend on the basis. In the basis (30), they 
are given in Appendix A. 



3.3. Algebraic equations for the density matrix 

Transformation properties of the coefficients a and A with 
respect to the reflection p z — > — p z or p' z —p' z allow us 
to split the system fl4(i|) into two independent subsystems 
for T) = 1 and n = — 1. Using the basis (30), we have 



and 



^Pr)ns,—s ^P—rj,ns,—S' 



(54) 



Thus it is sufficient to solve one of the subsystems: 

- ^ (A nSl s"Pns's2 + A* ns , S2 ip nSlS i) — 
rj'n' s's" 



(55) 



where 7 = 77' for s' = s", and 7=1 otherwise. In Eq. (55) 
we have set r\ = 1 in a, A, and y>. 

Retaining diagonal elements of ip and setting si = S2 
in Eq. (55) we recover the algebraic equations derived by 
Yakovlev (1984) in the distribution function formalism. 
However thus reduced system is not covariant with re- 
spect to the basis transformations ([32|) (note that the co- 
efficients a and A have been obtained from the matrix 
elements (52), and they undergo the transformations to- 
gether with the density matrix elements). The lack of co- 
variance leads to the dependence of the electron transport 
coefficients on the basis. For instance, Yakovlev (1984) ob- 
tained the difference up to 20% using different basis sets 
for some particular electron gas parameters. Our density 
matrix formalism makes the complete system ( 40 ) covari- 
ant with respect to the basis transformations. The trace 
of the density matrix is invariant, therefore changing the 
basis does no more affect $(e) (Eq. (42)) and the kinetic 
coefficients. 

3.4- Applicability range 

Let us discuss briefly the validity of the assumptions which 
led us to the algebraic system (55), for the neutron star 
crust conditions. First, we have assumed diagonality of the 
density matrix in p x and p z . This property is practically 
exact. Indeed, the delta- functions 8(q\ x — q2x) Hfliz — q2z) 
in Eq. (51) enter the right-hand side of Eq. ( |47| ) through 
Eq. (49) . Since p K i K > is diagonal in p' x and p' z , the right- 
hand side of Eq. ( jffl ) virtually contains S(pi x — P2x) 5{p\ z — 
P2z)- Thus the diagonality of the density matrix in p x 
and p z is not affected by collisions. Note that the delta- 
functions in Eq. (51) appeared due to the infinitely large 
volume assumed in averaging over ro . If the volume were 
finite, a non-diagonality in the momenta occurred in the 
band Ap ~ h/L, in agreement with the uncertainty prin- 
ciple. An actual value of L is restricted by the condition of 
spatial uniformity of the considered bulk of matter. The 
broadening Ap is negligible compared with the typical 
momentum difference between adjacent Landau orbitals, 
provided that L » h/(m c cVb), which condition is always 
satisfied in the magnetized neutron star crust. 
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Another assumption is concerned with the diagonality 
in the Landau numbers. The last term in Eq. (49) shows 
that this property is not conserved in collisions. This effect 
can be estimated with the aid of the uncertainty principle. 
The non-diagonality should occur within the collisional 
band Ae ~ Ti/r, where r is the effective relaxation time. 
The density matrix is diagonal in n if Ae < fwj* B . The 
electron cyclotron energy Tilob ranges from 10 2 to 10 4 Ry 
(where Ry =13.6 eV is the Rydberg energy), for magnetic 
fields from 10 11 to 10 13 G. For estimating Ae, we can use 
the non-magnetic relaxation times (27), (28) (cf. Yakovlev 
1984). If Z ~ 26, T ~ 10 7 K, and the electrons are 
mildly relativistic, we obtain %/t ~ Ae ~ 30 Ry < Tuob- 
This estimate relates to the collisional broadening. Other 
broadening mechanisms, not treated here, are due to non- 
elasticity of scattering and deviations from the Born ap- 
proximation. Yakovlev (1984) argued that the two latter 
types of broadening seem to be unimportant under the 
considered physical conditions. 

While considering the spin number s = ±1, one may 
notice that the spin-polarization density matrix p V ns's is 
reduced compared to the general case of 4 x 4 matrix 
(Bcrcstctskij ct al. 1982). It is because we restrict the basis 
to the electron bispinors Xns, thus neglecting an admix- 
ture of positron states. Of course, this restriction is well 
justified at the non-relativistic temperatures of interest. 

Finally, the assumed degeneracy in the spin number is 
not exact, owing to the quantum-electrodynamical cor- 
rections. The corrections split each Landau level into 
two sublevels (e.g., Landau & Lifshitz 1982). According 
to the Schwinger formula, the splitting energy is about 
(a/2w)hu} B ~ (0.1 - 10) Ry < Ae, for B = 10 11 - 10 13 G. 
Therefore the collisional width of the sublevels exceeds 
their separation, and they can be considered as degener- 
ate. 

Thus Eq. (55) is valid for the typical conditions in neu- 
tron star crusts. 

3.5. Numerical examples and discussion 

We have computed the first-order correction (p nslS2 e to 
the density matrix and its trace $(e) for a wide range 
of magnetic fields from 10 11 to 10 14 G. The coefficients 
in Eq. (55) have been calculated using the formulae from 
Appendices A and B. The system (55) has been solved 
using the LU factorization code by Fletcher (1988). We 
have calculated the density matrix and also the electron 
distribution function in two representations, with = 
and cf> = (j) n (Eqs. (30), @, and (§|)). The coefficients of 
the algebraic systems for the distribution functions have 
been taken from Yakovlev (1984), and the integrals Qi 
which enter these coefficients have been calculated using 
the formulae of Appendix B. 

At relatively low field strengths, the density matrix 
formalism gives practically the same results as the distri- 
bution function approach with the basis (30) (0 = 0), in 



agreement with anticipation of Yakovlev (1984). The in- 
accuracy of the latter approach is about 0.1% for B ^ 
10 12 G. The distribution function results based on the al- 
ternative "fixed helicity" representation (0 = <p n ), em- 
ployed by Hernquist (1984) and Schaaf (1988), deviate 
from the exact results by several per cent. When n e = 1, 
the deviations reach 3% for the phonon scattering and 8% 
for the Coulomb scattering. The deviations decrease at 
higher n e . 

The discrepancies between different representations 
are larger for stronger fields. Figure [l] shows the results 
for the iron crust (Z = 26) with magnetic fields B = 10 13 
and 10 14 G. The upper panels display the function $(e) 
for the electron energies sufficient to occupy up to ten 
Landau levels. One can see strong quantum oscillations 
of ^(e). Sharp dips at integer values of v are caused by 
singularities of the electron density of states due to the 
magnetic quantization, when electrons start to populate 
new Landau levels. 

The lower panels show relative errors of the distri- 
bution function approximation for two basis choices dis- 
cussed above. The discrepancies reach up to 20% at B = 
10 14 G for the Coulomb scattering, in agreement with re- 
sults of Yakovlev (1984). However, the energy dependence 
of the relative errors is non-monotonous, and a simple in- 
terpolation between the approximate results proposed by 
Yakovlev (1984) is not very accurate. The basis (30) pro- 
vides higher accuracy in most cases, but the alternative 
basis ([32]), (||^) becomes more appropriate at relativistic 
energies for phonon scattering. 

Generally, the discrepancies between different rep- 
resentations decrease with increasing energy. We have 
checked numerically, that relative errors do not exceed 1% 
for v > 50 even at B = 10 14 G. Calculations for such high 
v are facilitated by semiclassical formulae for the integrals 
Qi derived in Appendix B. 

Our calculations show that one should use the density 
matrix formalism for energies 1 < v <, 5 and relativistic 
magnetic fields b 1 , if one needs to keep an error within a 
few percent. Otherwise the distribution function approach 
with the basis (30) provides the desired accuracy. 

For a better comparison with the results of Yakovlev 
(1984) and Hernquist (1984), the examples in Fig.[l] have 
been calculated with the scattering potentials used by 
these authors. In particular, the screening length in the 
liquid regime has been determined by Eqs. (7) and (8). In 
the solid regime, we have adopted the conditions when 
the Debye - Waller factor can be neglected. Under these 
assumptions, $(e) does not depend on temperature. The 
effect of finite temperatures on ^(e) is generally small in 
the liquid regime. For the solid regime, the effect is dis- 
cussed in the next section. 
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4. High-temperature conductivity enhancement in 
the solid phase 

As seen from Eqs. (^) - (|TT|), the Debye - Waller reduc- 
tion of electron-phonon scattering becomes important at 
sufficiently high temperatures. Its influence on the elec- 
tron transport in dense non-magnetized stellar matter has 
been considered by Itoh et al. (1984a, 1993). The effect 
is most pronounced near the melting point. In this case 
it can be easily estimated. Let us put p — p$, where 
Pa = ?i(37r 2 n ) 1 / 3 is the field-free Fermi momentum. Tak- 
ing into account Eq. ([To|), we can write w in Eq. ( ^6|) as 
w = (4/3)u_ 2 (97rZ/4)273/r, which gives w = 560/r for 
u_2 = 13, Z = 26. Near the melting point T ~ 180, 
and w ~ 3. Then the reduction factors are i?i ~ 0.3 and 
i?2 ~ 0.2. According to Eq. (28), this leads to the largest 
increase of the relaxation time by a factor of 3.3 for non- 
relativistic electrons and 2.5 for ultrarelativistic electrons, 
in qualitative agreement with numerical results of Itoh et 
al. (1984a, 1993). 

However the Debye - Waller factor has been neglected 
so far while studying the electron transport in quantizing 
magnetic fields (Yakovlev 1984, Hernquist 1984, Schaaf 
1988, Van Riper 1988). In this case, contrary to the non- 
magnetic one, the longitudinal momentum p z is quan- 
tized into the discrete values mcqP^n (see Sect. fO]), and 
only scattering events with changing p z do contribute 
to the transport processes. Therefore there exists a low- 
est collisional momentum transfer Ap m ; n (e). According 
to Eqs. (All) and (A12) of Appendix A, the scattering 
rate acquires an additional exponential reduction argu- 
ment C£min = (fTAp m i n /?i) 2 /3. The smallest momentum 
difference reaches Ap m [ n ~ 2\f2m c huB just below the first 
Landau threshold, and Ap m j n ~ m c huiB/po just before 
each next threshold. Therefore the exponent argument is 
C£min ~ w for n e = 0, and C£ min ~ w(m c haj B /pl) 2 for high 
n e . In the latter case, the additional exponent argument 
becomes small and unimportant. However, it is significant 
at low n e . 

Figure ^ shows the Debye - Waller enhancement of 
the electron relaxation time in quantizing magnetic fields. 
For each T, B, and electron energy e we have put [i = e 
and calculated the electron number density from Eq. (||) 
to determine the Debye - Waller exponent argument. The 
fitting formulae presented in Appendix C facilitate this 
calculation. Then the coefficients of the algebraic system 
(55) were computed, using the formulae of Appendices 
A, B, and the system was solved to obtain the function 
$(e), which is proportional to the effective relaxation time 
r(e), Eq. (42). Its ratio to the same function calculated 
without the Debye - Waller factor is plotted against the 
Landau variable v in Fig. |^. This ratio reveals quantum 
oscillations in phase with oscillations of the function <&(e) 
itself, thus increasing the amplitude of the latter ones. For 
comparison, the non-magnetic enhancement factor (i?i — 
(i? 2 /2)(w/c) 2 ) -1 (cf. Eq. (28)) is shown by dashes. As seen 



from Fig. g, the quantizing field makes the Debye - Waller 
enhancement much stronger. When the electrons populate 
one or two Landau levels, the relaxation time is enhanced 
by a factor of up to 30 (compared with the factor of about 
3 in the non-magnetic case). 

5. Summary 

We have presented the theory of transport properties of 
degenerate electrons along quantizing magnetic fields in 
neutron star crusts. Our results are advanced, compared 
to the previously known ones, in three respects. 

First, a kinetic equation for the spin polarization den- 
sity matrix of electrons is derived. The solution of this 
equation provides a justification of the standard approach 
based on the kinetic equation for the electron distribution 
function. The present results are compared with two ver- 
sions of the standard approach used previously by different 
authors. 

For non-relativistic magnetic fields, B & 10 13 G, our 
results confirm the arguments of Yakovlev (1984) that the 
standard approach which employs basic functions with 
fixed spin z-projection is the most appropriate in the non- 
relativistic limit. The fixed-helicity basic functions used 
by other authors lead to small inaccuracies which however 
seem to be insignificant in astrophysical implications. 

For stronger fields, B <; 10 13 G, the inaccuracies of the 
traditional approach increase up to 20% when density is 
rather low and the electrons occupy low-lying Landau lev- 
els. The density-matrix results allow us to choose the most 
appropriate version of the standard approach. If density is 
higher and the electrons populate many Landau levels, the 
difference between various approaches becomes negligible. 

Secondly, we have taken into account the cumulative 
effect of the Debye - Waller factor with the magnetic quan- 
tization. In the non-magnetic case, this factor can increase 
the thermal and electric conductivities by a factor of 3 just 
below the melting temperature. We show that the mag- 
netic quantization can enhance the effect by an order of 
magnitude. 

Thirdly, we have derived semiclassical expressions for 
some intermediate integrals which enter the system of 
equations either for the density matrix or for a distribu- 
tion function. These expressions provide fast and accurate 
calculation of the relaxation time for large number of oc- 
cupied Landau levels. 

In this paper we have calculated the kernel function $ 
which should undergo further thermal averaging, Eq. ([l9]) , 
to determine the longitudinal electric and thermal conduc- 
tivities and thermopower. We shall consider this averaging 
and astrophysical implications of the developed theory in 
the subsequent paper. 
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Appendix A: coefficients in equations for the den- 
sity matrix 

While deriving Eq. (53), let us first reduce Eq. ( ^ ) to Eq. (^). 
For this purpose, let us sum Eq^lTf over p% x , p 2z , n 2 and p% z 
taking into account Eqs. @, (B6), (49) - (52). This yields 



(If 



rj' n' p' z s' s" 



s" Pr)ns" S2 (P»)+ 



i 



E 

rj' n' p' z s' s" 



(-^Cs' s i,s' s" Pr/ns" S2 i.P z ) 



2tt 

n 



^2 J]] Ks'si,s"s 2 Pri'n's's"(j>z) S(ui), 

rj' n' p' z s' s" 



(Al) 



where w = (e' — e)/h and 



Si ,8" S- 



\p' z ;n,p z ) = 



"i 



(2tt) 2 L 2 



q*dq y x 



\U (ci)\ 2 M*, 31 (n, p' z ; n,p z ; u)M a „ S2 (n ,p' z \n,p z ;u). (A2) 

It is convenient to introduce new integration variable u = (qcji + 
q 2 ,) (Jm/2 and the dimensionless function 



v(u,£) 



(hcam) 1 \U(q)\ 



V27T, 

For the Coulomb scattering, 

v ( U ,0 = ( U + ?)~\ £ = {Qz + rs 2 ) dm/2, 



(A3) 



(A4) 



where r s is the screening length, and Tiq z — p' z — p z . For the 
scattering on phonons, 



v{u,Q = (u + g) 1/2 e w , £ = q z a 2 a /2. 



(A5) 



The sum over p' z in Eq. (Al) can be converted into the sum 
over r\ and the integral over e . The delta- function 5{u) elimi- 
nates the integration in the first and the last terms of Eq. (Al). 
Finally we arrive at Eq. (53) with 



,(E,E), 



R,e -Arjnsi S2 / ^"r}nsiS2,f] f n f s f s' ; 

r} 1 n f s' 



(A6) 
(A7) 



Im An 



ry n s 



v 



E , 



dE' 
E' — E 



"rjns \ S2 ,r}' n' , 



,(E,E'), 



(A8) 



where E n / = vT+ 2bn' 



^rjnsi S2 n 1 s' s" {_E j E ) 



EE' 
\PP'\ 



du v (u, £) x 



Ms> 31 {n ,p' z ;n,p z ; u)M a „ S2 (n ,p' z ;n,p z ;u), 



(A9) 



and P = P V n(E) is defined in Sect. 3.1 According to Eqs. (52) 



and (30), the matrix elements M a / a contain the integrals 



i dy 

Ti-n'iy/am - qxa m /2)Ti n (y/a m + q x a m /2) e lgyV — = 

dm 

[(n')\/nl\ 1/2 e- u/2 u (n - n ' )/2 L n J n '(u) =/„„,(«), (A10) 



where I nn /(u) is a Laguerre function (Sokolov & Ternov 1968). 
Therefore M s i s can be expressed in terms of the functions 



h(£,n',n,m) 



Inn' ( u )In-l,n' -l( u ) _f( u -|-f) 



(All) 
du, (A12) 



and Qi = Q 2 (£,n' — l,n — l,m). Here m = 1 and £ = 
2ry/(3a 2 n ) for the scattering on phonons; m = 2 and £ = 
for the Coulomb scattering. The functions Qi, Q 2 , and Q 3 
have been introduced by Yakovlev (1984) for the particular- 
case of C, = 0. The properties of these functions are analysed 
in Appendix B. 

Using the basis (30), we obtain 



^nnsiS2,f]' n' s' s" \E , E ) 



S(s u S 2 \s's" 



4(E + 1)(E> + 1)\PP'\' 



(A13) 



where 



5(1, 1|1, 1) = g Qi + 4b 2 nnQ 2 + AbgV¥i~' Q 3 , (A14) 
5(-l, — 1 1 — 1, -1) = g 2 Q 2 + Atfnn'Qx + 4bgVmi/Qs, (A15) 
5(1,1|-1,-1)= 

2&(n'P 2 Qi + nP' 2 Q 2 -2\WPP'Q 3 ), (A16) 
5(-l,-l|l,l)= 

2b{n'P 2 Q 2 +nP ,2 Q 1 -2Vrm>PP'Q 3 ), (A17) 
S(-l, 1|1, 1) = 5(1, -1|1, 1) = -gV^P'Qx + 

2bn'V2bnPQ 2 + V2bn'(gP - 2bnP')Q 3 , (A18) 
5(1, 1| - 1, 1) = 5(1, 1|1, -1) = -gV2~briPQ 1 + 

2bnV2lm'P'Q 2 + V2lm(gP' -2bnP)Q 3 , (A19) 
5(-l, 1| - 1, -1) = 5(1, -1| - 1, -1) = gV2b^P'Q 2 - 

2bn'V2bnPQ 1 - V2bn'{gP - 2bnP')Q 3 , (A20) 
5(-l, -1| - 1, 1) = 5(-l, —1)1, -1) = gV2briPQ 2 - 

2bnV2bn'P'Q 1 - V2tm(gP' - 2bnP)Q 3 , (A21) 
5(-l,l|-l,l) = 5(l,-l|l,-l) = 

2bgV¥i 7 (Q 1 + Q 2 ) + (4b 2 nn' + g 2 )Q 3 , (A22) 
5(-l,l|l,-l) = 5(1,-1|-1,1) = 

2fe\WPP'(Qi + Q 2 ) - 2b (n'P 2 + nP' 2 )Q 3 , (A23) 

and g=(E + l) (E' + 1) + PP'. 

The first two pairs of coefficients (Eqs. (A14), (A15), and 
(A16), (A17)), are related to diagonal density matrix elements 
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for electron transitions without and with spin flip, respectively, where F is a Gauss hypergeometric function. Then we obtain 
At E — E' they reproduce equations of Yakovlev (1984). The 
remaining equations (A18) - (A23) present the coefficients at 

off-diagonal elements of the correction tp to the density matrix Q2+j((,, ri , n, m) = S ■ 
in the system (55). 



1/2 



[nl(n'y.(n-j)\(n'-j){ 
k\ (n-j- k)\ (ri -j - k)\ 



Appendix B: integrals Qi(£, n', n, m) 

For an efficient computation of the coefficients (A6) - (A8) in 
Eqs. (55), let us consider the properties of the functions Qi 
(Eqs. (All), (A12)). Since Qi(£,ri ,n,m) = Qi(£,n,ri , m), we 
assume ri — n > without loss of generality. 

Using the polynomial representation for the Laguerre poly- 
nomials, the integrals Qi can be expressed as 



Q 2 (£,ri ,n,m) 



2n 

E 

3=0 



(-1)^ 



k 



Cn' nk^-n' n ,j — k " 



(ri -n+ j)\Q 2 {E,,n' -n + j,0,m), 

2n-l 

Qa{£,,n' ,n,m) = } (-1) 3 } c n > nk c n ,„j_ k —= x 



j=0 k 

(ri - n + j)\ <2 2 (£,n' —n + j,0,m), 

(1 + 0~"'e £ £„, +1 + 
e« (E n ,(0-E n , +1 (0), 



Oa(€,n',0,l) 
Qa(£,n',0,2) 



(Bl) 



(B2) 
(B3) 
(B4) 



where Ei(£) is an integral exponent which is easily calculated 
(Abramowitz & Stegun 1972), 



y/ n'\n\ 



k\(n-k)\(ri -n + k)\ 



(B5) 



are the coefficients of the polynomial expansion of a Laguerre 
function I n i n (u), and summation index k runs from max(0, j — 
ri) to min(n, j). However this method fails for n <; 15 due to 
the exponentially increasing round-off errors (Hernquist 1984, 
Schaaf 1988). Direct integration (Hernquist 1984) is not very 
efficient because of rapid oscillations of the Laguerre functions 
at large n. The direct integration is especially undesirable for 
the density matrix computation, since an outer integration over 
E' is required in Eq. (A8). 

To avoid the above difficulties, we propose to use the trans- 
formation 



(to-1)! 



-(« + ?) x 



dx 



(B6) 



in Eqs. (All), (A12), and change then the integration or- 
der. Let us use the equalities (Gradshtein & Ryzhik 1965, 
Abramowitz and Stegun 1972) 



Jo 



u e >L n (w)L„_j (u)du 



(ri + n~jy. x 2n - j 



(ri - n)\ (n - j)\n\ (1 + x) n '+ n -'+ 1 
F y—n + j, — n; — ri — n + j; 1 — 1/a; 2 ) 



(B7) 



x ] ^ (ri)\(ri - j)\ 

(l + x) n + n '-j +1 2-*> (n-j-k)\(ri-j -k)\(j + k)\ k\ ' 

k — 



2k+j+m-l 



(m-l)\(J + k)\ J (1 +x) n '+ n ~i+ 1 



dx. 



(B8) 



The main advantage of the representation (B8) is that all terms 
are positive monotonous functions of £. 

Using Eq. (B8), we can easily derive various asymptotes. 
In particular, if £ = and £ — > oo, then from the asymptotic 
properties of the integral in Eq. (B8) we obtain 



Q2+j(£,n',n, m) 



(n-j)\(n'-j)\ 



1/2 



(j + m-l)\ 
j\(m-l)\ 



(l/?r +J [1 - (n + ri - j + 1) (j + m)/i + . 



(B9) 



The functions Qi are finite for ri 7^ n, and they diverge at 
£ — > for ri = n: 



An 

M£,n,n,l) = -ln£- 7 - V — + 

m, 



(ri) 2 v-^ (2n - 2m)! (2m - 1)! 



E 



(2n)\ ^ [(n-m)!m!] 2 

m=l 



+ O(0, 



(BIO) 



Q 3 (0 n, n, 1) = - In £ - 7 - V" h 



n!(n-l)!^ (2ra - 2m - 1)! (2m - 1)! 



E 



(2n - 1)! ^ (n-m)\(n-m - 1)! (m!) 2 

m — 1 

2n + l 



+ O(0. (BH) 



Q 2 (e,n J n,2)=r 1 + (2n + l) 

E 



(ri.) 2 sr^ (2n - 2m + 1)! (2m - 2)! 
(2n)\ 2^ [(n-m)!m!] 2 



in^+y -+■ 

^ m 

m = 2 

+ 0(0, 



+ 



(B12) 



Q 3 (Z,n,n,2)=C 1 +2n 



Zn 

in?+y 1+. 



+ 



ri. (n - 1)! x ^ (2n-2m)!(2m-2)! 
(2n - 1)! 2-^ (n-m)\(n-m- 1)! (m!) 2 



+ O(0- (B13) 



However, collecting the terms with ri = n in Eq. (55) , we see 
that only the finite terms of Eqs. (B10)-(B13) do contribute 
to the coefficients of the algebraic system, while contributions 
from the divergent terms (£ _1 and ln£) are mutually compen- 
sated in the real parts of the coefficients (Eqs. (A6), (A7)), as 
well as in the integrals (A8) for the imaginary parts. Therefore 
the coefficients of the density matrix equations are essentially 
finite, except for the singularities at the Landau threshol ds du e 
to the factor \PP'\ in the denominators of Eqs. (A9), (A1S). 
These remaining singularities stem from divergencies of the 
electron density of states at the Landau thresholds, and they 
are responsible for the qua ntum oscillations of the relaxation 
time discussed in Sect. 3.5. 
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Finally, consider the case n > 1. Let us make use of the 
semiclassical approximation of a Laguerre function averaged 
over oscillations (Kaminker & Yakovlev 1981): 



where e n = e n (0) is denned by Eq. (Fy) with p z — 0, and 



(B14) 



where 111,2 = (v« T V™) 2 - Using again the transformation 
(B6), we express the integrals Qi as 



1 r cxp (-ICyfnFn costf) 
Q 2 (Z,n',n,l) = ± / — — ; ' di9 (B15) 



(B16) 



71 Jo n ' + n + f + 2v n'ra cos 1? 

2a;\/ n'n 



dxe ^io 



(n' + n + 0y fK+n+f) '"V^' + n + C 
4Vn'nQ 3 ((;,n ,n, 1) = (n + n') [Q 2 (£,n', w , 1) + 
Q 2 (£, n' - 1, n - 1, 1)] + £ [Q 2 (£, n', n - 1, 1)+ 
Q 2 (£, n' - 1, n, 1)] - exp (~C(n' + rc + £ - 1)) x 

Jo (2C v / «'(^-l)) + /o (2C V(n' - 1) 
n', n, 2) = -{d/d® Qi(£, ri, n, 1), 



(B17) 
(B18) 



where Jo is a modified Bessel function. In deriving Eq. (B17), 
we have used the relationship (Kaminker & Yakovlev 1981) 

AVrvnJ I nn >{u) I n _ 1}n ,_ 1 (u) = (n+ ri) (-Z" 2 n '( u ) + 

/n-l,n'-l(«)) -« (-^n-l,n'(") +-fn,n'-l(w)) ■ (B19) 

The integral (B16) can be easily calculated using polynomial 
approximations for Io(x) (Abramowitz & Stegun 1972). At 
2V n'n/(n' + n + f) < 0.95, however, direct integration of 
Eq. (B15) is more efficient, since only a few integration mesh 
points are needed then. If £ = 0, the semiclassical approxima- 
tion allows us to express Qi in elementary functions: 



Q 2 (e,n',n,l) = [( U1 + (u2+£)n 1/2 , 

Q 2 (£, n', n, 2) = (£ + n + n') + f) (u 2 + ^}~ 3/2 , 

4Vnn' Qi(^, ri, n, 1) = 

(n + ri) (Q 2 (£, n', n, 1) + Q 2 (£, n' - 1, n - 1, 1)) + 
€ (Qa(£, n', n - 1, 1) + Q 2 fc, n' - 1, n, 1)) - 2, 

4v / ^7Q 3 (£,n',n,2) = -((n + n')£ + (n -n) 2 ) x 
(<?!(£, n' - 1, n, 1) + Ql(e, n', n - 1, 1)) + 
(n + ri) (Q 2 (£, n', n, 2) + Q 2 (^, n' - 1, n - 1, 2)). 



(B20) 
(B21) 



(B22) 



(B23) 



We have checked numerically that errors in calculating $(e) 
do not exceed 3%, if we substitute the semiclassical formulae 
(B15) - (B17) and (B20) - (B23) into Eqs. (A14) - (A23) for 
n > 10. The semiclassical approximation makes the computa- 
tion much faster. 

Appendix C: a fitting formula to electron density 
Equation (|^) can be represented as 



2(-Kh) 2 c 



(CI) 



Evidently, 



dn e m, 
~d/7 ~ 2(rrn) 



k B T 'k^Tj ' 



(C2) 



(C3) 



where F'(x,y) = dF{x,y) / dx. 

In the non-relativistic limit we have y ^> 1 and 1 « j 
(or x < 0). Then F(x,y) reduces to y/y/2J r _ 1 / 2 (x), where 
J--i/ 2 {x) is a Fermi integral, for which some useful approxi- 
mations have been presented by Antia (1993). In the opposite 
case of y <C 1 we have F(x, y) — In (1 + e x ). 

In the general case of arbitrary y, the following approxima- 
tion is proposed: 



F(x,y) = \n(l + e x 



1 + y + g + c(y) p(g) VCT2^ 
l + f + c( V )6tt) 



(C4) 



where £ = ln(l + exp[a; — xo(y)]), 
x {y) = {l + c iy ^)~\ 
c(y)=c 3 y Ci , 

a(0 = A /?/2 + (ai + a 2 £ 2 )Vf J 
6(C) = 1 + a 3 V? + + a 2 £ 3 , 

with the numerical parameters ci = 0.623, c 2 = 1.6031, C3 = 
0.9422, c 4 = 1.7262, 01 = 0.103, a 2 = 0.043, a 3 = 0.0802, and 
a 4 = 0.2944. At y > 1, the right-hand side of Eq. (04) simpli- 
fies to ^ V£ + 2y a(f)/6(£), with £ = ln(l + e 3: ), and depends es- 
sentially only on the four constants a% . An error of the approx- 
imation (C4) reaches maximum of 0.57% at x = 2.9, y = 5.9. 
We have checked additionally, that the ^-derivative of this ap- 
proximation fits the function F'(x,y), which enters Eq. (p3j), 
with a maximum relative error of 2% . 
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Figure captions 



Fig. 1. Function $(e) given by Eq. (42). Landau variable v is defined by Eq. (^). 
Top panels: <E>(e) in the density matrix approach. 

Bottom panels: Relative errors A<J>/$ = ("l > appr/ c I > — 1) of the distribution function approximations $ appr (dotted line: basis 
(30); dot-and-dashed line: basis (^), (|3^)). Dashes show non- magnetic approximation, according to Eqs. (27), (28), and (41). 



Fig. 2. Ratios of the effective relaxation times TDw(e) to Tnon-Dw(e) calculated, respectively, with and without the Debye - 
Waller factor. 

The temperature T = T 7 x 10 7 K, with T 7 = 0.2, 0.5, 1, 2, 5, and 10. Solid lines - numerical results (Eqs. (41), (42)), dashed 
lines - non-magnetic case (Eqs. (27), (28)). The curves are plotted starting from the melting points. 
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